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Abstract 

Several isotropic, homogeneous cosmological models containing a 
self-interacting minimally coupled scalar field, a perfect fluid source 
and cosmological constant are solved. New exact, asymptotically sta- 
ble solutions with an inflationary regime or a final Priedmann stage 
are found for some simple, interesting potentials. It is shown that the 
fluid and the curvature may determine how these models evolve for 
large times. 



PACS 04.20 Jb, 98.80 Cq 



1 Introduction 



A self-interacting scalar field has been introduced in cosmological models 
as a matter source to the Einstein equations because, when dominated by 
the potential energy, it violates the strong energy condition and drives the 
universe into an inflationary period fl] , |§ , @] . Scalar fields also play a fun- 
damental role in unified theories of the strong, weak and electromagnetic 
interactions through the mechanism of spontaneous symmetry breaking. Of 
particular interest for cosmology are those symmetries which spontaneously 
broken today were restored in the hot early universe. Thus, several kinds 
of topological defects may have been produced 0,0. Further, a scalar field 
has been proposed as a kind of dark matter ||. 

Currently, there is no underlying principle that uniquely specifies the 
potential for the scalar field and many proposals have been considered. Some 
were based in new particle physics and gravitational theories Q . Others were 
postulated ad-hoc to obtain the desired evolution. For instance, Ellis et al. [|7| 
have proposed a scheme to find the potential function for a given expansion 
of a Robertson- Walker universe. Also, a formalism has been proposed to 
reconstruct the potential from knowledge of tensor gravitational spectrum or 
the scalar density fluctuation spectrum [EJ, H. 

Notwithstanding the number of papers devoted to understand the dy- 
namics of the self-interacting scalar field in curved spacetimes, very little is 
known yet about exact solutions of these cosmological models. Most of these 
solutions correspond to models in a spatially flat Robertson- Walker metric 
which have no other source besides the scalar field [|TJ, |TD[, JTT[| , [O, JTJ], 



[PH ) |15| ) [[U|, |17|], [pi , fl9|| . Few exact solutions are known with spatial 



curvature PD| , 0, pifl , fl2"2fl , or a perfect fluid 0, p3| , p^[ ; and we are aware 
of none with both terms. The problem arises in the non-linearity of the sys- 
tem of differential equations for the scalar and gravitational fields. Usually, 
in studies of inflationary dynamics, the evolution is divided in several steps, 
and during each of them some kind of approximation is assumed to simplify 
the system of equations ||, 0. 

In section 2 of this paper we show a procedure to reduce to quadra- 
tures the Einstein-scalar field equations in a Robertson- Walker metric for 
an arbitrary potential, a perfect baryotropic fluid source and a cosmological 
constant. This allows us to address the issue of whether the scalar field is al- 
ways the dominating force in driving the evolution of the universe or the fluid 
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and curvature terms may also play a role. We analyse in section 3 general 
properties of the solutions and we study their asymptotic stability by means 
of the method of Lyapunov |25[] . We consider also their structural stability. 
Our procedure is illustrated in section 4 with several examples. First, we use 
it to rederive and generalize a few solutions for the well known exponential 
potential, and then we present some new exact solutions for several simple, 
interesting potentials. Some of them are suitable for cosmological models 
based on the chaotic inflationary scenario. No a priori assumption like a 
slowly varying field is required to perform the calculations, so that we can 
check the validity of this assumption. Finally, the conclusions are stated in 
section 5. 



2 The Model 

We wish to investigate the evolution of a universe filed with a scalar field and 
a perfect fluid. The scalar field <fi has a self-interaction potential V(4>) and is 
minimally coupled to gravity, so that it obeys the Klein-Gordon equation 



, dV 
□0+ — = 



(1) 



The perfect fluid has four-velocity and its pressure p and energy density 
p are related by the equation of state 



p = ( 7 - l)p (2) 

with a constant adiabatic index < 7 < 2. Thus, we must solve equation 
(lip together with the Einstein equations 



Rik - ^gikR + Ago* = T? k + T/ k 



(3) 



We are using units such that c = 811G = 1, A is the cosmological constant 
and 



T; 



ik 



:k — 9ik 



T; 



ik 



(p + p)uiU k - pg ik 



(4) 
(5) 



are the stress-energy tensors of the field and the fluid. 
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In a Robertson- Walker metric 



ds 2 



dt 2 - a 2 (t) 



dr 2 



1 — kr 2 



+ r 2 (d6 2 + sin 2 9d4> 2 ) 



(6) 



with scale factor a(t) and curvature parameter k = 0, ±1, equations (^) and 
(H) become 



+ 3H(j) + -y- = 



3iJ 2 = ^-0 2 + V{4>) + p — 3— + A 
2 cr 



(7) 



(8) 



where the dot means d/dt, H = a/ a and <fi = <p(t). Also, from the conserva- 
tion of (H), p = Po/ a? 1 where po > is a constant. 

To integrate this system we make the change of variables dt = a 3 drj in 



d 2 <f) s dV () 
drf d(j) 



and we write 



V[cj>{a)\ = 

Thus we obtain the first integral of (|7]) 



F(a) 



(9) 



(10) 



where C is an arbitrary integration constant. Then, using the scale factor 
as the independent variable and equations (|8|) and flTT|), we have reduced the 
problem to quadratures: 



At = V%J 

a<p = Vg[- 

J a 



da 
a 



6 r F C & . 

da— + — +p-3— + A 



a" J a a 6 a 2 



-1/2 



—F + 6 J daF/a + C 
6 / daF/a + C + pa 6 - 3ka A + Aa 6 



1/2 



(12) 
(13) 
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where At = t — to, A0 = — O and to, 0o are two other arbitrary integration 
constants. 

Equations ([TI]), ( JT2] ) and ([□]) allows us to obtain a solution a(t), 4>(t) for 
the potential V(4>), given a function F(a), the parameters p , 7; k and A, 
and a choice of the primitive / daF/a and the integration constant C. 



General Behavior of Solutions 



In general the set of solutions that arise from ( Jig ) (|13| ) contains complex 
solutions without physical sense that must be rejected. So, to obtain physical 
solutions we require that the functions 

6 f F C k 
G{a) = — da- + — + p - 3— + A (14) 



a J a a nZ 



F 6 f F C 
L ( a ) = --*+-J da ~ + -K ( 15 ) 



a 6 a v J a a v 

are nonnegative on an interval (a\, 02), < a\ < ci2 < 00. Since G(a) = 3H 2 
and L(a) = <f) 2 /2, a sufficient condition is that a > 0, a 7^ and 0^0. 
To relax these conditions we need to study the behavior of the scale factor, 
which is determined by the function G(a). Thus, a solution a(t) is monotonic 
(d 7^ 0) when G(a) > but it is bounded if G(ao) = for < a < 00. A 
growing monotonic solution starts from a = and expands without bound. 
A solution is singular when a(t) vanishes at a finite time, that is when the 
integral 

m , x f a2 da . . 

T( ai ,a 2 )= ^_ (16) 



converges in the limit a% — *■ 0. Let us assume that in this limit the leading 
behaviors of F(a) and / daF/a are Ka n and Ka n /n, where K > and n/0 
are constants. So G(a) ~ K/a a where if > is another constant and a is the 
maximum of the exponents in (0), namely (6 — n, 6, 37, 2). Then solutions 
are singular when a > 0, which occurs when any of the following conditions 
are met: C > 0, po > 0, k = — 1 or n < 6. Singular solutions have particle 
horizons when a > 2, that is when C>0orn<4or7> 2/3. Similarly, 
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assuming F(a) ~ a n for a — > oo, we find that no solution diverges at a finite 
time. 

For bounded solutions, we observe that the scale factor reaches a at a 
finite time to when T (0,1,0,2) converges in the limit ai — > ao- Let us assume 
that G(a) ~ .^(a 5 — Oq)^, /3 > and 5^0 constants, in a neighborhood 
of Oq. If /5 < 2 then oq is a local extremum of a(t) and the evolution is 
time-symmetric about t . Otherwise, for (3 > 2, a(i) — > a when |£| — > 00. 

Provided that the requirements for a physical solution are met, the free- 
dom to select F(a) reflects the freedom to choose the potential. However, 
only in restricted cases the functions a(t), <p(t) and V(4>) arise in closed ex- 
plicit form. 

Though the solution ( JT2] ) flT3| ) of the system of equations (|7|) (Q) has three 
constants of integration, it does not provide the general solution for a given 
potential. Rather it provides a special solution for each member of a family 
of potentials. So, for each potential we need to study the stability of this 
solution. 

For open or flat models such that V(4>) has a local minimum at <p m and 
A e ff = V(4) m )+A > 0, we can study the stability of solutions with asymptotic 
behavior <p(t) — > <fi m . First we note that <fi = <p m and a(t) given by 

Ai ^/?('-! +A -P < 17 > 

is a classical solution of (|7|) and (|8|) with the effective cosmological A e ff > 
[p6[| (a(t) may not be given by ( |P7D when A e ff = 0). 

To study the asymptotic stability of a solution it becomes convenient 
to turn it into a fixed point. For the models under consideration a(t) is 
monotonic we may use a as the independent variable. Then we can restrict 
the analysis to the phase space (4>, (ft') and study the stability of the point 
(4>mi 0)- We find that the energy density of the field 

P4> = v(<j>) + l -4> 2 (18) 

is a suitable a Lyapunov function as it satisfies p<j,((j), , a) > —A and p$ = 
—3aH 2 <f)' 2 < in a neighborhood of (<f) m ,0). So this point is an attractor 
and any solution such that (f> — > (f> m for a — > 00 (equivalently t — > 00) is 
asymptotically stable. In this limit, when A e ff > 0, equation (|T7| ) gives the 
leading behavior of a(t). Since 
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"=-¥ 2 -b" + ^ < 19) 

is negative we note that these solutions have a deflationary behavior. 

On the other hand, in closed models some solutions are non-monotonic. 
However, for any solution such that the scale factor is monotonic starting 
from a given time, we may apply the same analysis as before about its stabil- 
ity. Also, for k = 1 we have the scalar field analog of the Einstein Universe 
p6fl . In effect using (0), (||) and ( |I9"D we find for a nonconstant potential a 
unique static solution a = a E > and (p = <f) E , provided that V'^e) = 0, 

V(<p E ) = 3 -^-A (20) 

l a E 

and 7po — 2a^ r ~ 2 . In the case 7 = 2/3 there is no particular equilibrium 
density and there is a static solution for any a E - To study the stability of 
the static solution we consider small perturbations about it. In the linear 
regime the system decouples and we find that the Einstein solution is stable 
when 7 < 2/3 and V"(4>e) > 0, and it is unstable otherwise. 

Recently it has been considered the question whether small variations in 
the form of the potential induce small variations in the scale factor and vice 
versa. A model with this property is called structurally stable or "rigid" 
2% . In the case that k = A = p = 0, the rigidity of inflationary models has 
been shown ||28|| . We wish to investigate whether this property is also valid 
without these restrictions. Let us consider a family of functions F(a) labeled 
by a parameter r. If F(a,r) is a continuous function and a > 0, using (|T0|) 
and (0), a small change 5r induces a small change SV in V{<p). Also, using 



(|T2"D, it induces a small change 5a in a(t). Clearly, small neighbourhoods in 



the space of scale factors and in the space of potentials map continuously 
into each other. Thus we conclude that models where F is continuous are 
structurally stable. 



4 Examples 

In this section we obtain several new exact solutions for simple interesting 
potentials and we concentrate mainly on models with spacial curvature and 
a fluid source since very few exact solutions are known for them. 



7 



4.1 Exponential and exponential- like potentials 

We consider first potentials of the form 



V{4>) = J Bexp(-(xA0) 



(21) 



Exponential potentials like (^Tj) have been extensively considered by many 
authors, mainly in the context of "power-law" inflationary models. We red- 
erive and generalize exact solutions for this potential. For this we take 

F(a) = Ba s (22) 

where B and s are constants, we choose / daFj a = Ba s / s with C = and we 
obtain a real solution when < s < 6 for B > 0. Since s = 6 corresponds to 
a constant potential, we consider this case no further. For s < 6, the energy 
density of the scalar field redshifts as a perfect fluid with equation of state 
P<j> = (l<p - l)pf P<t> = Po^/a™, where p 04 , = 6B/s and n = 3^ = 6 - s [|]. 
Using fllT]) and ([12]) we find a power-law solution of the form 



a(t) = (KAt) x A(f) = Mlna 



(23) 



with some constants K, M and A which are functions of the parameters of 
the system. 

Case 1. p Q = k = A = 0. 



a 



(6 



a/2 



K 



M = (6 - s 



.1/2 



This solution has been found by Lucchin et al. ||10|| , its asymptotic stability 
has been analysed by Halliwell P0| and the general solution for this model 



has been found by Salopek et al. ||15 
Case 2. p = A = 0, s = 4, k 



0. 



a 



B 



[B-2k) 



1 or k 

1/2 



A 



M 



1 with B > 2. 

B_ xV2 
2 

1/2 



K = 

2B 
B-2k 
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This solution has been found by Ellis et al. [0 (see also [5D| 



for k = 



-I)- 



Case 3. k = A = 0, s = Q — 37. 



a 



(6 



1 + 



Pqs 
6B 



1/2 



K 



6B 



3s 



s Po \ 1/2 Q-s 



A 



6 



2 

37 



M 



65(6 



nl/2 



6fi 



This solution has been found by Ratra et. al || in the case 7 = 1, and 
we extend it here to arbitrary 7. It generalizes case 1 for a fluid source 
and for 7 > 2/3 it is not inflationary and has particle horizons. It is usually 
assumed that an exponential potential with a small enough leads to a power- 
law inflationary evolution (A > 1), so that any matter term becomes rapidly 
dominated by the scalar field energy density. However in our exact solution, 
which is asymptotically stable, the energy densities of the field and the fluid 
keep a constant ratio along the evolution. 

Another interesting generalization of case 1 arises from a nonvanishing 
constant C, which yields a family of exponential-like potentials. First we 
note that C > for a real solution. Then, using (|2~2"D , fll3l) and ( |10|) we find 
the potential in parametric form 



A4>(a) = — (w + 2, 

s { V 6 



6 - s\ 1/2 



In 



a s + Cl ) 1/2 + (a s + c 2 



M 2 



2 In 
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(a s + Cl ) 1/2 + (a s + c 2 ) 1/2 



(24) 



V{a) = Ba s - 6 (25) 

where C\ = sC/((6 — s)B) and c 2 = sC/(6B). The leading behavior of this 
potential is ( pip with a = — (6 — s) 1 / 2 for <\> — > 00, while a = — (6 — n)/v6 
for — > —00. 

When C > we find that 



At 



/ s \V2 /1 3 3 a s 



27 1 - 1 7 
5 5 Ci 



(26) 
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where %Fx is the hypergeometric function [2j|. These solutions also begin at 
a singularity, this time like a ~ At 1 / 3 for At — > 0, while their asymptotic 
behavior is like case 1 for t — > oo. Though these solutions have some initial 
conditions fixed, we note that they are asymptotically stable. 

A common framework to solve equations (^) , (H) in discussions of inflation 
is the " slow-roll" approximation. However, as it has not been required by our 
procedure of integration, we can use it to investigate the limitations imposed 
by this approximation. Following Copeland et al. [§], and using equations 
(ED (0) we calculate the slow-roll parameters 

-- ^ -l- = (27) 



V + 02/2 C + 6 J daF/a V 6 J a s + c 2 



4> 12F - aF' - 36 / daF/a - 6C _ (s - 6) a s - 6ci 
V = Jl4>~ 2(C + 6/daF/a-F) 2 (a s + ci) ^ ^ 

The parameter e measures the relative contribution of the field's kinetic 
energy to its total energy density and rj measures the ratio of the field's 
acceleration relative to the friction term. The slow-roll approximation is 
valid when | e |<C 1 and | rj |<C 1. These requirements impose constrains on 
the form of the potential and the value of the initial conditions. 

In this case, for C — 0, the slow- roll constrain reads 6 — s 1, that is, 
this approximation is valid only when the potential has a very mild slope and 
A > 1. Moreover, for C > 0, the evolution near the singularity is dominated 
by the kinetic energy. Thus, in general, these solutions cannot be obtained 
under the slow-roll approximation. 

4.2 Cosh and cosh-like potentials 

We consider now potentials of the form 

V{(j)) = V [cosh (aA(j))} q (29) 

with Vq > 0. This potential has a minimum at the origin when q > or a 
maximum when q < with a value Vq. For large \<j)\ it grows or vanishes 
exponentially. The case q < appears when reconstructing the potential in 
the assumption of small amplitude of the tensor gravitational wave spectrum 
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M. For q > we have the stable solution = 0o and a(t) given by ( |TT| ) with 
A e // = Vo. We consider bellow other solutions for this potential. 

Case 1. F(a) = Ba s (b + a s ) n , where B > 0, b > 0, s and n are constants 
such that s{n + 1) = 6, s ^ 0, s ^ 6, and p = k = A = C = 0. 



V = B a 



We obtain the time evolution for arbitrary s in implicit form: 

< 3 '3 3 . 3 



At 



a {6 6 6 a 



Acf) 



2VQ 



arccoth ( 1 + 



s \ 1/2 



(30) 



(31) 



It is suitable for further analysis to consider explicit expressions for the lead- 
ing behaviour of a(t) and <j)(t) in the limit of small and large a s /b 



a(t) ~ (VWb^'At) 1 



/3 



a(t) ~ exp 

, 2 v 1/2 



B\ 1/2 A ' 
— At 



3 

In At, 



a s < 6 



(32) 
(33) 

A0~ (J) In At, a s <6 (34) 

A0^-^exp[-j(fy /2 Atl, a s »6 (35) 

(36) 

^-Htt^J (37) 

shows that the slow-roll approximation is valid only for a s ^> b if \s\ <C 1. So 
these are new solutions which cannot be obtained under such an approxima- 
tion. 



The calculation of the slow-roll parameters 

b 



b + a s 
1 / 66 + sa s 
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We may describe the behavior of these solutions as follows. When < 
s < 6, a(t) begins at a singularity with the behavior (p2|) for At — > and 
is asymptotically de Sitter like ( |3"3"D for At — > oo. The evolution of the field 
<j)(t) starts dominated by a diverging kinetic energy and settles down at the 
minimum of the potential for large times driving the exponential expansion. 
When s < the evolution begins like (^) in the remote past and then 
deflates into the asymptotic Friedmann behavior ([32]) in the far future. The 
field begins rolling down very slowly from the maximum of the potential. 
Then it gains kinetic energy that decays much more slowly than the potential 
energy. As it becomes dominant, it drives the deflationary era and makes the 
field grow without bound. We provide in this way an exact realization of the 
deflationary universe scenario proposed recently by Spokoiny J3"0| . We note 
also that for either sign of s, these solutions are stable. Finally, when s > 6, 
the scale factor behaves the same as the case < s < 6, but the field begins 
its evolution dominated by the kinetic energy and approaches very slowly to 
the maximum of the potential for large times. This is an unstable solution. 



Case 2. F(a) = Ba s , A = C = 0, s = 6 - 37, 0<7<2, k = 1. 

1 / 2B 



a 




3 V2-7 



+ Po 



37 
2-37 



2 I 67 
The time evolution is given by 

2 /w\V2 



2 + f (2-7) 



1/2 



67 



37-2 



1 37 3 



1 



At = — I - a 37/2 2 Fi . . 

37 \b) V 2 67-4' 2 2-37 



1 / o _o\!/2 

arctanh ( 1 — 00a 7 



uoa 



37-2 



(38) 
(39) 



where 
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For 0<7<2/3(g<0), the scale factor has a bounce and its asymptotic 
behavior for large times is t 2 ^ 3/y \ For 2/3<7<2(g>0), the evolution 
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starts from a singularity with the leading behavior At 2 ^ 3 ^ , it reaches a max- 
imum and then collapses again. We note that for the physically interesting 
case of radiation 7 = 4/3, the solution takes a simple closed form: 



a(t) 



B+ P -± 
3 



At 



1/2 



(40) 



Case 3. F(a) = Ba s , Po = A = C = 0, < s < 6, s ^ 4, k = l. 

4-s 2(6 - 1 



6=8 

4-s 



The time evolution is given by 



a 



2^6 



q 



4-s 



KB) "6^7 2 1 l2'2(s-4)'2(s-4)'TB~ 



1 , / sa 4 " s 

arctanh 1 

a I 25 



1/2 



(41) 



(42) 



For 4 < s < 6 (q < 0), the scale factor has a bounce and its asymptotic 
behavior for large times is £ 2 /( 6_s ). For < s < 4 (q > 0), the scale factor 
starts from a singularity with the leading behavior At 2 ^ 6 ~ s \ it reaches a 
maximum and then collapses again. 

We show here another potential that have qualitatively the form (|29|): 



V{<f>) = B 



(2+37) 

24 



where e = ±1 and 



eexp(crA0) 



3e|exp(cxA0) 



^ exp(2crA</>) 



37-4 
2-37 



(43) 



2-37 

>A 4 - 3 ^) 



(44) 



They arise in the case that F(a) = Ba s , A = C = 0, s = 2 + 37 and 
k — 1. There is a critical energy density coefficient p c = 3(2 + 37)/(85) and 
for po < p c we find a cosh-like potential with c/ > when either 7 < 2/3 
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(e = 1), 2/3 < 7 < 4/3 (e = -1) or 2/3 < 7 < 4/3 (e = 1), this time 
q < 0. The scale factor starts from a singularity with the leading behavior 
At 2/(4-3 7 \ it reaches a maximum and then recollapses again. 

We see from our examples that we cannot discard a priori the influence 
of a curvature term in the determination of the evolution, since a potential 
that leads in the flat case to an asymptotically expanding scale factor yields 
recollapsing solutions when k — 1. 



4.3 Sinh and sinh-like potentials 

We consider potentials of the form 

V((f)) = V \smh(aA(j))\ q (45) 

with Vq > 0. This potential has a vanishing minimum at the origin when 
q > or diverges when q < 0. For large \<f>\ it grows or vanishes exponentially. 
This kind of potentials satisfy, for suitable values of q > 0, the requirements 
usually imposed in the chaotic inflationary scenario M. We study several 
exact solutions for this potential. 

Case 1. F(a) = Ba s , A = C = 0, s = 6 - 37, < 7 < 2, k = -1. 

'1 / IB 



a 



V = B 
37 - 2 f 1 



3 \2-7 



+ Po 



37 

2-37 



2 [67 
The time evolution is given by 



At = a 2 F 1 - 



2 + §(2- 7 ) 



1/2 



67 



37 



,2-3^ 



1 + 



2' 2 -37' 2-37' u 

1/2 



Acf) = arccoth f 1 + 

a v 



una 



37-2 



(46) 



(47) 



The evolution of the scale factor for these models starts from a singularity 
as At Xl for At — » and grows monotonically with asymptotic behavior t X2 
for t — > 00. In this case, for < 7 < 2/3, Ai = 1 and A2 = 2/(37) > 1) while 
for 2/3 < 7 < 2, Ai = 2/(37) < 1 and A 2 = 1. 



14 



Case 2. F(a) = Ba s , k = A = C = 0, < s < 6, 6 - 37 - s ^ 0. 



V = B 



spo\ 



-3-f-s 



a = 



6 — 37 — s 



.6BJ 2x/6^ 
The time evolution is given by 



2(6- s) 
6 — 37 — s 



At= (2s) 1/2 ^^^^ (1 6-s 3(6 -27-3) _Po£ /7 6_ 37 _; 



Vfi7 6-s 2Fl V2' 2(6 -37- s)' 2(6-37 -sj' 65° 



A0 = -larccoth f 1 + ^a 6 " 3 ^^ 

(7 \ DB / 



1/2 



(48) 
(49) 



For 6 - 37 - s > 0, q > 0, Ai = 2/(6 - s) and A 2 = 2/(37). Then, Ai < A 2 
but this is not an inflationary solution unless 7 < 2/3. For 6 — 37 — s < 0, 
q < 0, Ai = 2/(37) an d ^2 = 2/(6 — s). So Ai < A 2 , and this model may 
describe the passage from a primordial radiation-dominated era to a power- 
law inflationary stage. 

Case 3. F(a) = Ba s , A = C = 0, s = 4, k = 0, -1 or k = 1 if B > 2. 



Po V2 



2 

37-2 



V = £ 

The time evolution is given by 

n 

At 



37-2/I k " 
2 I 2 ~ 5 



1/2 



9 = 



2-37 



(f-fc) 1 "" 'V2'2-37' 2-37 
A0 = -^arccoth (l + cua 2 " 37 ) 1/2 



(50) 



(51) 



For 7 < 2/3 g > and we find Ai = 1 and A 2 = 2/(37) > 1- For 7 > 2/3 
q < 0, Ai = 2/(37) and A 2 = 1. 

Case 4. F(a) = Ba s , Po = A = C = 0, < s < 6, s ^ 4, Jfe = -1. 



25/ 



(i-s 
4-s 



er = 



4-s 
2 V / 6^B 



9 = 



2(6 -s) 
4-s 
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The time evolution is given by 



A</> = - -arccoth (l + 7^ 4 " s ) " (53) 

For < s < 4 q > 0, A x = 2/(6 - s) < 1 and A 2 = 1. For 4 < s < 6, q < 0, 
Ai = 1 and A 2 = 2/(6 — s) > 1, so that the evolution is inflationary. 



The expression fl43|) (e = 1) also yields potentials of the sinh form, with 
g > for po > Pc- They occur when = 1 and 7 < 2/3, or = — 1 and 
< 7 < 4/3, 7 7^ 2/3. The scale factor starts from a singularity and grows 
monotonically. For < 7 < 2/3 Ai = 2/(4 - 37) < 1 and A 2 = 2/(37) > 1, 
while for 2/3 < 7 < 4/3 Ai = 2/(37) < 1 and A 2 = 2/(4 - 37) > 1. In either 
case there are particle horizons and the evolution is inflationary. 

We note that all the solutions presented for potential (|4"5| ) are asymptot- 
ically stable. Many of them show the influence of the fluid source on their 
behavior, either near the singularity or for large times. 



5 Conclusions 

We have shown that the Einstein equations with a self-interacting minimally 
coupled scalar field, a perfect fluid source and cosmological constant can be 
reduced to quadratures in the Robert son- Walker metric . For that, the scale 
factor is taken as the independent variable and the potential is expressed in 
terms of it. This is fully justified as long as a(t) and <f>(t) are monotonic 
functions. For this class of cosmological models we have reobtained in a 
unified manner many exact solutions presented by other authors, and found 
several new exact solutions for some simple interesting potentials. We have 
investigated their behavior in detail, and shown that most of the solutions 
that do not recollapse in a finite time are asymptotically and structurally 
stable. 

For a given potential, the influence of a fluid source or the spatial curva- 
ture may dominate the evolution of the early universe. In effect, an asymp- 
totically inflationary behavior cannot occur when the energy density of the 
scalar field decays faster than the energy density of a fluid with 7 > 2/3. We 
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have shown asymptotically stable solutions such that both densities keep a 
constant ratio, or one of them dominates for large times. All the same, a neg- 
ative spatial curvature may lead to a linearly expanding scale factor. On the 
other hand, a positive curvature has frequently as a consequence bouncing 
solutions or solutions that undergo recollapse to a second singularity. In such 
cases, the energy density of the scalar field may increase, instead of rolling 
down the potential. 

We have studied the limitations imposed by the slow-roll approximation 
and most of our solutions cannot be obtained by means of this approximation. 
Calculations of the primordial spectrum of perturbations, based on recent 
measurements of the cosmological background radiation show that tensorial 
perturbations may have played an important role in the formation of cosmic 
structures. However, the amount of gravitational perturbations predicted by 
means of the slow- roll approximation is very small BT| . Thus exact solutions 



which lay outside the slow-roll regime may lead to an improved understanding 
of the evolution of the early universe. 
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